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Abstract 

The improved next-to-next-to- leading order (NNLO) QCD analysis of the exper- 
imental data of the CCFR collaboration for the xF% structure function is made. 
Theoretical ambiguities of the NNLO fits are estimated by means of the Pade re- 
summation technique, which was applied both in the expanded and non-expanded 
forms. The NNLO and the new N 3 LO a s (Q 2 ) MS-matching conditions are used. 
In the process of the fits we are taking into account the target mass corrections and 
the twist-4 l/Q 2 -terms. Our NNLO results for a s (Mz) values, extracted from the 
CCFR xF 3 data, are a s {M z ) = 0.118 ± 0M2(stat) ± 0.005(syst) ± 0.003(theory) pro- 
vided the twist-4 contributions are fixed through the infrared renormalon model and 
a s (M z ) = 0Al6±° o Q o ° Q 6 7 (stat) ± 0.005(syst) ± Q.00A(theory) provided the twist-4 terms 
are considered as the free parameters. It is shown that the extracted at the NNLO 
order x-shape of the twist-4 correction is almost unchanged after application of the 
Pade resummation approximant. 
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1. Introduction. 

Deep-inelastic lepton-nucleon scattering (DIS) belongs to the classical and continu- 
ously studying processes in the modern particle physics. The traditionally measurable 
characteristics of vN DIS are the SFs F2 and XF3. It should be stressed that the 
program of getting the information about the behavior of the SFs of vN DIS is among 
the aims of the experimental program of Fermilab Tevatron and CCFR/NuTeV collab- 
oration in particular. The CCFRR collaboration started to study the vN scattering 
process over 1980 year §. The data for the SFs of vN DIS, obtained by the follower of 
the CCFRR collaboration, namely CCFR group, was distributed among the potential 
users in the beginning of 1997 while the final results of the original CCFR DGLAP 
H NLO analysis of this data was presented in the journal publication of Ref.Q. 

This experimental information was already used in the process of different NLO 
analysis, performed by CITEQ, MRST and GRV groups (see Refs.@, | § corre- 
spondingly). The subsequent steps of performing NLO and first NNLO analysis of the 
CCFR data with the help of the Jacobi polynomial - Mellin moments version of the 
DGLAP method were made in Refs. j7t- [ fl2[ (the definite stages in the development of 
this formalism are described in Refs.|R|^15fl) 

In the process of the analysis of Refs. p[|-||12|| the authors used the information about 
the NNLO corrections to the coefficient functions |l6| and the available analytical 
expressions for the NNLO corrections to the anomalous dimensions of the NS moments 
with n — 2, 4, 6, 8, 10 [jl7|, supplemented with the given in Ref. |1 n = 3, 5, 7, 9 similar 
numbers, obtained using the smooth interpolation procedure, which was previously 
proposed in Ref. |l8|] . Moreover, the attempts to extract the shape of the twist-4 
contributions and study the predictive abilities of the IR-renormalon (IRR) model of 
Ref. |l{| were performed (for the definite details of modeling the effects of the power- 
suppressed contributions to measurable physical quantities using the IRR language 
see Ref. 20 and the review of Ref. [^lj ) . 

However, the important question of the estimation of theoretical uncertainties of 
the NNLO analysis of the CCFR data of Ref. [fllf was still non-analyzed in detail. 
These uncertainties are determined by 

1) the differences in the definitions of a s (Q 2 ) matching conditions (see e.g. |22|, 
24|) which are responsible for penetrating into the energy region, characteristic for 

/ = 5 numbers of flavours, where the pole of the Z°-boson is manifesting itself; 

2) the incorporation into the condition of Ref.|22[ recalculated NNLO QCD cor- 
rections 1 25 1 and the newly calculated N 3 LO corrections, namely the 4- loop coefficient 
of the QCD /^-function |§6) and the N 3 LO-term |2^] in the matching condition of 
Ref.||; 

3) the consideration of theoretical uncertainties due to other non- calculated N 3 LO 
contributions into the coefficient functions and the anomalous dimensions function. 

This work is devoted to the analysis of the important problems outlined above and 
to the more detailed extraction of the values of a s (Mz) and the x-sh&pe of the twist- 
4 power-suppressed term at available orders of perturbative QCD with taking into 
account the effects enumerated above. We are supplementing the NNLO fits of Ref. |nj 
by the N 3 LO analysis, which is based on the application of the Pade resummation 
technique (for the review see Ref. |28[] ) , developed in QCD in the definite form in 
Refs.|2^, |3TJ| and considered previously as the possible method of fixing theoretical 
uncertainties in the analysis of DIS data in Ref. pl[| . It should be stressed that a 
posteriori this technique gives the results similar to those, obtained with the the help 
of the application of different methods of fixing scale-scheme dependence ambiguities 
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(compare the results of Ref. |32| with the results of Refs.[g9| |30| obtained using the 
Pade resummation technique). Thus, our analysis could be considered as the attempt 
to estimate perturbative QCD uncertainties beyond the NNLO level. Moreover, it 
could give us the hint whether the outcomes of the NNLO fits, related to perturbative 
and non-perturbative sectors, stay stable after the inclusion of the explicitly calculated 
and estimated N 3 LO QCD corrections. 

2. The theoretical background of the QCD analysis. 

Let us define the Mcllin moments for the NS SF xF 3 (x, Q 2 ): 



x'^F^x^^dx 



(1) 



where n = 2, 3, 4, .... The theoretical expression for these moments obey the following 
renormalization group equation 



9 nl A N 9 



r NS 



where A s — a s /(47r). The renormalization group functions are defined as 



^=(3(A S 

O/J, 



i>0 



/'- 



dlnZ, 



NS 



dfi 



(2) 



(3) 



where Z^ s are the renormalization constants of the corresponding NS operators. The 
solution of the renormalization group equation can be presented in the following form 



exp 



As(Q 2 ) An) 



As(Q 2 ) 



7ns( x ) 
0(x) 



dx 



(4) 



where M^ s (Qq) is the phenomenological quantity related to the factorization scale 
dependent factor. It can be parametrized through the parton distributions at fixed 
momentum transfer Qq as 



M» S {Ql)= / x^A{Ql)x b ^\l-xr^\l + 1 {Ql)x)dx 



(5) 



with 7 ^ or 7 = 0. In principle, following the models of parton distributions, used 
in Refs.||, |(J, one can add in the used model for the SF the term, proportional to y/x. 
However, since this term is important in the region of rather small x, we will neglect 
it in this our analysis. 

At the N 3 LO the expression for the coefficient function cl?c can be presented as 



C$l(A.) = 1 + C (1 \n)A s + C (2 \n)A 2 s + C^(n)A 3 s , 
while the corresponding expansion of the anomalous dimensions term is 



exp 



0(x) 



-dx 



(A 5 (<2 2 )) 7jvs(n)/2/3 ° x AD(n s ) 



(0) 



(7) 
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where 

AD(n s ) = [l+p(n)A s (Q 2 ) + q(n)A s (Q 2 ) 2 + r(n)A s (Q 2 ) 3 ] (8) 
and p{n), q(n) and r(n) have the following form: 

P[n) ~2\ X O )(3 



(9) 

*(») - I (W) 2 + ^ + 7^W^P - 7&(")§) (10) 



r(n) = ± ( p{nf + 6p(n) q (n) + ^J^L , , , , 

falttlin) , /3?7$» $7#k») &7$» , Wifojfiliny 



0o @o Po A) Aj 

The coupling constant A S (Q 2 ) can be expressed in terms of the inverse powers of 
L = ln(Q 2 /Af^) as A? LO = A L ° + AA^ LO , A? NLO = A? LO + AA? NLO and 

Af LO = A^ NLO + AAf LO , where 

A L ° = (12) 
aa nlo = Piln{L) 



0*L 2 
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AA» NLO = ^ [/^n 2 (L) - #fn(L) + &A, - A?] (13) 

AAf LO = _Lj[$(_f„3(£) + \ln 2 {L) + 2fn(L) - \) (14) 

-3/3 /3i/3 2 Zn(L)+/3 2 ^] . 

Notice that in our normalization the numerical expressions for /3q, Ab 02 and /?3 
read 

/3 = 11 - 0.6667/ 

01 = 102 - 12.6667/ 

02 = 1428.50- 279.611/ + 6.01852/ 2 

3 = 29243.0 - 6946.30/ + 405.089/ 2 + 1.49931/ 3 (15) 



where the expression for 0^ was obtained in Ref.|26 . The inverse-log expansion for 
AAf LO , which incorporates the information about the coefficient 0a, was presented 
in Ref.|7|. 

Few words ought to be said about the used approximation for the anomalous 
dimension function jjfy(A s ). The analytical expression for its one- loop coefficient is 
well-known: 7^(n) = (8/3)[4£" =1 (l/j) - 2/n(n + 1) - 3]. In the cases of both F 2 
and xF$ SFs the numerical expressions for jjp s (n)-coefficients are given in Table 1. 
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n 


fire c (fl) 


tSuc j? (n) 




/iV5 v y I L-L/ ij 


Tlrc(?l)|[n/9l 


1 


2.5575 


o 


306 6810 


? 


? 


2 


71.3744 


71.2410 


612 0598 


5259 


5114 


3 


100.8013 


100.7819 


837.4264 


6959 


6900 


4 


120.1446 


120.1401 


1005.8235 


8421 


8414 


5 


134.9049 


134.9035 


1135.8235 


9563 


9562 


6 


147.0028 


147.0023 


1242.0056 


10493 


10482 


7 


157.3323 


157.3321 


1334.0017 


11310 


11280 


8 


166.3862 


166.3861 


1417.4506 


12077 


12012 


9 


174.4683 


174.4682 


1493.5205 


12784 


12706 


10 


181.7808 


181.7808 


1559.0048 


13370 


13271 


11 


188.4662 


188.4662 


? 


? 


? 


12 


194.6293 


194.6293 


? 


? 


? 


13 


200.3496 


200.3496 


? 


? 


? 


14 


205.6891 


205.6891 


? 


? 


? 



Table 1. The used numerical expressions for the NLO and NNLO coefficients of the 
anomalous dimensions of the moments of the NS SFs at / = 4 number of flavours and 
the N 3 LO Pade estimates. 

These results are normalized to the world with / = 4 numbers of active flavours. 

(2) 

In the same Table the numerical expressions for 7jVs( rl )i used in the process of the 
fits, are presented. In the cases of n = 2, 4, 6, 8, 10 they are following from the explicit 
calculations of 7)^5 F2 (n)-terms Jl?]], normalized to / — 4, while the n — 3,5,7,9 
numbers were modeled using the smooth interpolation procedure, originally proposed 
in Ref. Note in advance, that since 7^ F3 (n)-coefficients differ from (11)- 
terms, though by the small additional contributions (for discussions see Ref. H), it 
would be interesting to verify the precision of the model for 7^5 used in the 
process of our NNLO xF^ fits, by the explicit analytical calculations. 

Let us now describe the procedure of fixing other theoretical uncertainties. After 
the work of Ref. [^9| it became rather popular to model the effects of the higher order 
terms of perturbative series in QCD using the expanded Pade approximants. 

In the framework of this technique the values of the terms (n) and r(n) could 
be expressed as 



Pade [1/1] : C (3) (n) = [C {2) (n)] 2 / C {1) (n) (16) 

r{n) = q(n) 2 / p(n) (17) 

Pade [0/2] : C (3) (n) = 2C (1) (n)C (2) (n) - [C (1) (n)} 3 (18) 

r(n) = 2p(n)q(n) - [p(n)] 3 (19) 

The numerical expressions for p(n) and q{n), obtained from the results of Table 1 
and definitions of Eqs.(9)-(11), together with the values of the coefficients C^ln) and 
C^(n) (which are coming from the calculations of Ref. |l6| ) , are presented in Table 
2. 
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11 


p(n) 


q(n) 


K n )l[i/i] 


K")l[0/2] 


CW( n ) 


C7W(n) 


C W (»)ln/H 


CW(»)| r o/ai 


2 


1.646 


4.232 


10.829 


9.476 


-1.778 


-47.472 


-1267.643 


174.408 


3 


1.941 


4.774 


11.738 


11.218 


1.667 


-12.715 


97.004 


-47.013 


4 


2.050 


5.546 


15.003 


14.123 


4.867 


37.117 


283.085 


246.009 


5 


2.115 


6.134 


17.790 


16.486 


7.748 


95.408 


1174.834 


1013.328 


6 


2.165 


6.595 


20.087 


18.407 


10.351 


158.291 


2420.569 


2167.903 


7 


2.210 


7.039 


22.421 


20.318 


12.722 


223.898 


3940.284 


3637.790 


8 


2.252 


7.525 


25.138 


22.471 


14.900 


290.884 


5678.657 


5360.371 


9 


2.294 


8.018 


28.027 


24.715 


16.915 


358.587 


7601.720 


7291.305 


10 


2.334 


8.375 


30.049 


26.382 


18.791 


426.442 


9677.390 


9391.308 



Table 2. The expressions for the NLO and NNLO QCD contributions, used in 
our fits, and the N 3 LO Pade estimates. 



In the same Table we also give the estimates for r(n) and C^(n), obtained using 
the expanded [1/1] and [0/2] Pade approximants formulae of Eqs.(16)-(19). In the last 
two columns of Table 1 the estimates for the N 3 LO contributions to the anomalous 
dimension function j^l(A g ), obtained with the help of the expanded [1/1] and [0/2] 
Pade approximants are presented. One can see that the results of applications of [1/1] 
and [0/2] Pade approximants for 7]vs( n ) are almost identical to each other. 

Using the numbers, presented in Table 1, one can construct Pade motivated ex- 
pressions for r(n) by substituting the estimates for 7^(n) into Eq.(ll). It should 
be stressed, that the obtained by this way estimates for r(n) will qualitatively agree 
with the ones, presented in Table 2 within the "Pade world" only, namely only in 
the case of application in Eq.(ll) of the [1/1] or [0/2] Pade estimate for the four-loop 
coefficient of the QCD /3-function (3%. However, in the case of / = 4 the direct appli- 
cation of the [1/1] and [0/2] Pade approximants underestimates the calculated value 
of /?3 by over the factor 2.5 (/33 1 [1/1] = 3216.66...; /?3|[o/2] = 3058.38..). In view of this 
the application of Eq.(ll) with the Pade estimated values of 7jvs( r an d the explicit 
expression for /^-coefficient are giving estimates of r(n), drastically different from the 
ones presented in Table 2 (for example, for the case of application of [0/2] Pade esti- 
mates it gives r(2) w 18.83,. ...,r(10) « 55.70). It is already known that the accuracy 
of the estimates of the N 3 LO coefficient of the QCD /3-function can be improved by 
some additional fits of polynomial dependence of /% on the number of flavours / and 
applying the asymptotic Pade approximant (APAP) formula |33| . Therefore, it might 
be interesting to think about the possibility of putting bold guess Pade estimates of 
N 3 LO contributions to 7ats(^4 s ) (see Table 1) on more solid background. The anal- 
ogous steps were already done in Ref. [Q in the case of the analysis of the status of 
N 3 LO Pade estimates for the anomalous dimension function of quark mass and the 
agreement with the explicitly calculated at this level of QCD the results of Ref.j35| 

turned out to be reasonable. One can hope, that the application of the similar proce- 

(3) 

dure for the APAP estimates of ^j^g-terms and the substitution of the results obtained 
into Eq.(ll) (together with the explicit expression for the /?3-term) might improve the 
agreement with the results, presented in Table 2, which at this stage we consider as 
the most stable results for the N 3 LO fits. 

It should be also stressed that the uncertainties of the values of r(n) are not so 
important, since the results of our fits will be more sensitive to the form of the Pade 
approximations predictions for the contributions into the coefficient function (namely 
C( 3 )(n)-terms). 
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From the results presented in Table 2 one can conclude, that the theoretical series 
for C^g for large n (n > 4), which corresponds to the behavior of xF 3 (x 1 Q 2 ) SF 
in the intermediate and large a;-region, probably have sign constant structure with 
asymptotically increasing positive coefficients. Therefore, the applications of the ex- 
panded [1/1] and [0/2] Pade approximants are giving theoretical estimates for the 
terms (n) (n > 4) , which in both cases have the same positive sign and the same 
order of magnitude. 

However, in the cases of n — 2, 3 our intuition does not give us the idea what 
might be the sign and order of magnitude of the third term in perturbative series 
CP S (A„) = 1 - 1.784, - 47.47 A? and Cj$j(4) = 1 + 1.674 - 12.714?. Indeed, in 
these two cases the manipulations with [1/1] and [0/2] Pade approximants are giving 
drastically different estimates for the terms (n), which in the cases of n — 2, 3 differ 
both by sign and value (see Table 2). These facts, and the structure of the NNLO 
perturbative series for C NS (A S ) especially, might indicate that this series is not yet 
in the asymptotic regime. Another possibility is that its coefficients do not have the 
(+l) n n! growth, but poses some zigzag structure, which is manifesting itself in the 
cases of definite perturbative series of quantum field theory models (for discussions 
see e.g. Ref. |3(j]). This might give the additional theoretical uncertainties of modeling 
higher-order perturbative QCD predictions for F%(x,Q 2 ) in the region of small x. 

In view of the questionable asymptotic behavior of the NNLO series for the coef- 
ficient functions of NS moments with low n (n = 2, 3) we are also using the idea of 
Ref. |3lj] and consider the results of applications of non-expanded Pade approximants 
in the process of the analysis of the DIS data. 

Let us remind that the corresponding non-expanded [1/1] Pade approximants can 
be defined as 

^""'■'«=ttK: (20) 

14- r (n) A 

C&Wlp/H - 1)7 (21) 

where = ([p(n)} 2 - q(n)J/p{n), b[ n) = -q(n)/p(n) and c[ n) = (jC«(n)] 2 - 

C( 2 )(n))/C«(n), d[ n) = -C«(n)/CW(n). 

The explicit expressions for the non-expanded [0/2] Pade approximants read: 



&i nJ 4 + ^ J 4 



^•"'^ " i + 4-^+4-^ (23) 

where b[ n) = -p(n), b { 2 n) = p(n) 2 - q(n), d ( { l) = -C«(n) and 4" } = [C {1] {n)] 2 - 
C^ 2 \n). Since we consider the applications of both [1/1] and [0/2] Pade approximants 
as the attempts to model the behavior of the perturbative series for the NS Mellin 
moments beyond the NNLO level, we will use in Eqs.(20)-(23) the N 3 LO expression 
for the coupling constant 4, defined through Eqs.(12)-(14). It is worth to mention 
here, that quite recently the expanded and non-expanded Pade approximants were 
successfully used for the study of the N 3 LO approximation of the ground state energy 
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in quantum mechanics p7\ and of the behavior of the /3-function for the quartic Higgs 
coupling in the Standard Electroweak Model |Q . 

The next step is the reconstruction of the structure function xF 3 (x 1 Q 2 ) with taking 
into account both target mass corrections and twist-4 terms. The reconstructed SF 
can be expressed as: 

xF 3 N —(x,Q 2 ) = x a (l-xf £ e^(x)Y,c { ;\a,(3)Mj + 2, X F 3 (Q 2 ) (24) 

n=0 j=0 

, Kx) 
Q 2 

where ©°'^ are the Jacobi polynomials and a, (3 are their parameters, fixed by the 
condition of the requirement of the minimization of the error of the reconstruction 
of the SF. In order to take into account the target mass corrections the Nachtamnn 
moments 

M n>xF3 -> M™g{Q 2 ) = f -^F 3 (x,Q 2 ) l±^^ , (25) 

j x yn -t- z) 

can be used, where £ = 2x/(l + V), V = + 4-M i 2 t[ci a; 2 /(2 2 and M nuc i is the mass of 
a nucleon. However, to simplify the analysis it is convenient to expand equation ( p5|) 
into a series in powers of M 2 ucl /Q 2 ||{|. Taking into account the order 0(M^ uc i/Q i ) 
corrections, we get 



«W 2 ) = M^ F3 (Q 2 ) + n( ;+ 2 1} ^ < w (g a ) (26) 
, (n + 2)(n + l)nMl ci . 

We have checked that the influence of the order 0(M^ uc[ /Q 4 ) terms in Eq.(26) to 
the outcomes of the concrete fits is very small. Therefore, in what follows we will use 
only the first two terms in the r.h.s. of Eq.(26). 

The form of the twist-4 contributions h(x) in Eq.(19) was first fixed as 

N maj! (n) 

h{x)=x«{l-xf e ^WE c S n '( a ^) M /fwW 2 ) (27) 

n=0 j=0 

where Cj(a, 0) are the polynoms, which contain a and /3-dependent Euler T-functions 
and 

Ml R x %(Q 2 ) = C(n)M^ F3 (Q 2 )A' 2 + 0(±) (28) 

with A 2 taken as the free parameter and C(n) defined following the IRR model es- 
timates of Ref.(l| as C{n) = -n - 4 + 2/(n + 1) + 4/(n + 2) + 4Si(n) (Si(n) = 
X^=i l/i)' ^ should be stressed that the appearance of the multiplicative QCD 
expression M^^ F3 (Q 2 ) in Eq.(28), generally speaking different from the intrinsic co- 
efficients function of the twist-4 contribution, is leading to definite theoretical un- 
certainties in the contributions of higher-order QCD corrections to the twist-4 part 
of xF 3 (x,Q 2 ). This could provide the additional theoretical errors in the studies of 
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the the status of the IRR-model predictions for the twist-4 terms at the NNLO and 
beyond. 

In order to study this question at more definite theoretical level it is instructive to 
consider the function h(x) as the free parameters of the fits, not related to IRR-model 
estimates. 

We will estimate the uncertainties of the values of A-^L, a s (Mz) IRR-model pa- 
rameter A 2 and the twist-4 function h(x) due to the inclusion of the definite explicitly 
calculated N 3 LO QCD corrections in the fits and modeling other ones with the help of 
the Pade resummation technique. Our aim will be also the verification of the stability 
of the results of Ref . and the study of the possible reason of the peculiar behavior 
of h(x) for xF 3 SF, discovered in Ref.|ll| at the NNLO level. 

3 (a). The results of the analysis of the experimental data: the extrac- 
tion of AilL vs a s value. 

MS 

The definite results for our NLO and NNLO fits, made for the case of / = 4 
number of active flavours, are presented in Table 1 of Ref.|^], where the values of 
the parton distribution parameters A, 6, c, 7 7^0 are also given. In the present Section 
we study the influence of the Pade-motivated estimates of the N 3 LO expressions for 
anomalous dimensions and corresponding coefficient function (written down both in 
the expanded and non-expanded forms) to the results of the fits , which are resulting 
in the extractions of the parameter A.^L (and thus ot(Mz)), and of the common factor 

A 2 of the IRR model. 

It should be stressed that despite the general theoretical preference of applica- 
tions of the diagonal Pade approximants (for the recent analysis see e.g. Ref. po|| ) 
the N 3 LO [1/1] Pade approximant description of the CCFR'97 experimental data 
turned out to be not acceptable in our case, since it produces rather high value of 
\ 2 : \ 2 /nep > 2 (where nep = 86 is the number of the experimental points, taken 
into account in the case of the cut Q 2 > 5 GeV 2 ). However, the application of [0/2] 
Pade approximants produced reasonable results. We think that the nonapplicability 
of the [1/1] Pade method in the process of fitting CCFR xF 3 data using the Jacobi 
polynomial approach can be related to the manifestation of rather large value of the 
ratio [C (3) (2)/C (2) (2)]|[i/i] in the expression for NS moment 

The similar effect of the preference of the [0/2] Pade approximant analysis over the 
[1/1] one was found in Ref.|5(| in the case of the comparison of the QCD theoretical 
predictions for the polarized Bjorken sum rule (which are closely related to the QCD 
predictions for the first moment of the xF 3 SF, namely for the Gross-Llewellyn Smith 
sum rule) with the available experimental data. 

The results for A^L, obtained at the LO, NLO, NNLO and N 3 LO, modeled by 
the expanded and non-expanded Pade approximants, are presented in Table 3 in the 
cases of both 7 7^0 and 7 = 0. 

Looking carefully on Table 3 we arrive to the following conclusions: 

• Our fits demonstrate that the effects of the NNLO perturbative QCD contri- 
butions are important in the analysis of the CCFR data. Indeed, for different 
Q 2 -cuts they are diminishing the values of the QCD scale parameter A^- by the 
contribution, which is varying in the range over 50- 120 MeV, provided twist-4 
corrections are taken into account through the IRR model of Ref. jli| ; 

• The values of Ai-tL., which are coming from the fits with taking into account the 
[0/2] Pade estimates (both in the expanded and non-expanded variants) turn out 
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7 - free 


7 = 0- fixed 


Q 2 > 




A 2 (H1) 


X 2 /points 


A Ws i - McY > 


A 2 (tiL) (GeV )) 


X 2 /points 


o Uev 














T /■ \ 


266±37 


- 


113.2/86 


231±37 


- 


126.8/86 




436±56 


-0.33±0.06 


82.8/86 


367±59 


-0.30±0.06 


102.0/86 


ATT ( \ 

1NLU 


341±41 


- 


87.1/86 


278±33 


- 


109.3/86 




371±31 


-0.12±0.05 


81.8/86 


301±36 


-0.10±0.06 


105.7/86 


NATT n 

IN IN LU 


293±29 


— 


78.4/86 


258±28 


- 


100.6/86 




293±29 


-0.01±0.05 


78.4/86 


258±29 


-0.01±0.05 


100.5/86 


1\t3t r\ 
IN ±jU 


305±32 


- 


79.0/86 


306±33 


- 


88.5/86 


(n.e.) 


308±32 


-0.03±0.05 


78.5/86 


310±33 


-0.04±0.05 


87.8/86 


i\t3t r\ 


293±28 


- 


79.0/86 


295±30 


- 


87.1/86 




294±28 


-0.02±0.05 


78.8/86 


297±29 


-0.03±0.05 


86.8/86 


1U uev 














LU 


287±37 


- 


77.7/63 


279±40 


- 


79.6/63 




531±74 


-0.52±0.11 


58.0/63 


503±97 


-0.49±0.15 


61.7/63 


ATT \ 

1NLU 


350±36 




64.4/63 


328±39 




69.1/63 




/ion i c c 

439±55 


AO/) 1 A 1 A 

-U.24±U.10 


CO £ 1 GO 

58.6/63 


392±173 


A 1 A I A AT 

-0.19±0.27 


er O 1 GO 

65.3/63 


NNLO 


308±34 




58.7/63 


285±34 




68.1/63 




313±37 


-0.03±0.09 


58.6/63 


285±38 


-0.01±0.09 


68.1/63 


N 3 LO 


303±33 




60.4/63 


303±33 




63.0/63 


(n.e.) 


309±36 


-0.03±0.09 


60.3/63 


310±37 


-0.04±0.09 


62.8/63 


N 3 LO 


295±30 




59.9/63 


296±31 




62.2/63 




297±33 


-0.01±0.09 


59.8/63 


298±34 


-0.02±0.09 


62.2/63 


15 GeV' z 














LO 


319±41 




58.8/50 


318±41 




58.9/50 




531±57 


-0.57±0.13 


50.2/50 


517±112 


-0.54±0.27 


50.9/50 


NLO 


365±40 




53.0/50 


355±44 




54.1/50 




441±44 


-0.25±0.13 


50.9/50 


408±43 


-0.19±0.13 


52.8/50 


NNLO 


314±37 




50.9/50 


294±36 




56.7/50 




308±45 


0.03±0.14 


50.8/50 


279±44 


0.09±0.15 


56.3/50 


N 3 LO 


304±36 




52.8/50 


303±36 




54.3/50 


(n.e.) 


297±43 


0.05±0.14 


52.7/50 


294±43 


0.05±0.14 


54.1/50 


N 3 LO 


296±33 




52.3/50 


295±33 




53.8/50 




286±38 


0.07±0.14 


52.0/50 


283±39 


0.08±0.14 


53.4/50 



Table 3. The results of the extractions of the parameter A^L and the IRR coefficient A 2 , (in 
GeV 2 ) defined in Eq.(23), from LO, NLO, NNLO and N 3 LO non-expanded (n.e.) and expanded 
Pade fits of CCFR'97 data. 
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to be really nonsensitive to choosing the Q 2 -cut of the data, fixation of the value 
of 7 and thus incorporation of (1 + 7x)-factor in the parton distribution model. 
This in its turn can indicate that the change of the used by us model xF${x, Qq) = 
A(Q 2 )x»^(l-x) c ^(l+ 7 (Q 2 ))x) to xF 3 (x, Q 2 ) = A(Q 2 )x b ^(l+ 7 (Q 2 )x) + 
e(Ql)y/x), used in the MRST and GRV fits, should not affect significantly the 
results obtained; 

The large errors in the definite LO and NLO results for Ai^L, presented in Table 
3, are reflecting the correlations of these values with the errors of parton distri- 
butions parameters, which will be not considered in this paper and presented by 
us elsewhere ; 

For all Q 2 -cuts and values of the parameter 7 the applications of the N 3 LO 
fits performed with the help of the expanded [0/2] Pade approximants technique 
result in the smaller value of A^j and slightly smaller x 2 value than in the case of 
application of the non-expanded [0/2] Pade approximants. In our future studies 
we will consider the results of application of both expanded and non-expanded 
Pade approximations. 

For all Q 2 -cuts the expanded N 3 LO results for A^L is smaller then the similar 
results for non-expanded Pade fits by 10-15 M eV. This effect is coming from the 
expansion of Pade approximants in Taylor series. The difference in the results 
of the expanded and non-expanded Pade approximants can be considered as the 
estimate of the part of theoretical error of the N 3 LO results. 



Using the LO, NLO, NNLO and N 3 LO variants of the rigorous MS'-scheme match- 
ing conditions, derived in Ref.|2^] following the lines of Ref.[^2), we tr OTi - ■<'<•>■"■ \' ,! 
values through the threshold of the production of the fifth flavour M 5 



we transform Ai 

MS 

mi, (where 



rrib is the 6-quark pole mass) and obtain the related values of A^L with the help of 
the following equation: 



I3l +1 ln 



A 



2 

'MS 



(/3 /+1 - %)L h 



A 



(/) 2 

'MS 



+$NLO + $NNLO + Sn 3 



LO 



(29) 



Snlo 



Pi 



1% 



P{_ 

pi 



InLy 



Pi 



Pt 1 



In 



PT 

pi 



(30) 



'NNLO 



Piu 



+ 



pi 



f+i 



pi 
pi 



Pi ^Pi Pi ' 



Pi^ 

Pi +1 



pi ' 

— — ^2 



(31) 
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?N 3 LO 



(Plu, 



+ 



Pi 



f r 



pI 



pr 1 & 



pr \pi 



i 

" 2 
p{ 



Pi 



Pi 



4^ •/• 

pi 



(32) 



Pi 



Pt 1 



Pi 



/+1 



/+1 



/+1 



Pi 



f+1 



Pl +1 



PI 

Pi 

Pi 



Pl +1 



Pi 



/+1 



pi 


lnh\ 


3 pi + \ 


Pl\ 


pi +i 


Pi) 




-c 3 



where C2 = —7/24 was calculated in Ref.|2[| and the analytic expression for C3, 
namely C* 3 = -(80507/27648)C(3)-(2/3)C(2)((l/3)ln2+l)-58933/124416+(//9)[C(2)- 
2479/3456] was recently found in Rcf.|27|. Here /?/ are the coefficients of the 

/3-function with / (/ + 1) numbers of active flavours, Lh — /n(A/| +1 /A^L 2 ) and Mf+\ 
is the threshold of the production of the quark of (/ + l)-th flavour. In our analysis 
we should take / = 4 and m,b ~ 4.6 GeV. 

In the case of the non-zero values of the twist-4 function h(x) ^ the results of 
the LO, NLO, NNLO and N 3 LO fits, made both in the expanded and non-expanded 
versions of the [0/2] Pade motivated approach are presented in Table 4. 

It should be stressed that we are considering the outcomes of our N 3 LO ap- 
proximated fits as the rate of theoretical uncertainties of the NNLO results in the 
same manner like the results of the NNLO analysis will be considered as the mea- 
sure of theoretical uncertainties of the NLO results. In particular, we will intro- 
duce the characteristic deviations A NNLO = \(A^L) n3lo - (a£L) nnlo \, A nlo = 



(4) \NNLO 



lv MS J 



(A. 



(4) \NLO\ 
MS 



MS' 



It is worth to emphasize that the results for the Pade motivated N 3 LO results 
for A^L and thus a a (Mz) became closer to the NNLO ones (provided the statistical 
error bars are taken into account, see Table 1). Moreover, the difference A NNLO — 
\(A^-) n3lo - (A^§-) NNLO \, is drastically smaller then the NLO correction term 



(4) \NNLO 



lv MS' 



(K^L) NLO \ 

K MS I 



The similar tendency A 



NNLO 



« A 



NLO 



IS 



taking place in the case of the fits without twist-4 corrections. These observed prop- 
erties indicate the reduction of the theoretical errors due to cutting the analyzed 
perturbation series at the different orders. 

It is known that the inclusion of the higher-order perturbative QCD corrections 
into the comparison with the experimental data is decreasing the scale-scheme theoret- 
ical errors of the results for A^L and thus a s (M z ) (see e.g. Refs.Q fo|). Among 
the ways of probing the scale-scheme uncertainties are the scheme-invariant methods, 
namely the principle of minimal sensitivity, the effective charges approach, ( which 
is known to be identical to the scheme-invariant perturbation theory) and the BLM 
approach (for the review of these methods see e.g. Rcf. |^2j). The scheme- invariant 
methods were already used to estimate the effects of the unknown higher order cor- 
rections in SFs (see Ref.Q], where a strong decrease of the value of the QCD scale 
parameter was found in the process of the scheme-invariant fit of the experimental 
data for F 2 SF ) and to try to predict the unknown at present N 3 LO corrections to 
the definite physical quantities fl33], and DIS sum rules among others. Note that the 
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predictions of Ref.|52| are in agreement with the results of applications of the Pade 
resummation technique (see Ref.[^9j). Therefore, we can conclude that the applica- 
tion of the methods of the Pade approximants should lead to the reduction of the 
scale-scheme dependence uncertainties of the values of a s (Mz) in the analysis of the 
CCFR data. 

We are presenting now the values of a s (Mz), extracted from the fits of the CCFR'97 
experimental data for the xF^ SF, obtained with fixing twist-4 contribution through 
the IRR model of Rcf.g§: 

NLO HT of [19] a s (M z ) = 0.121 ± 0.002(stat) (33) 
±0.005(syst) ± 0.005(theory) 

NNLO HT of [19] a s (M z ) = 0.118 ± 0.002(stat) (34) 
±0.005{syst) ± 0.003 (theory) . 

The results of the extractions of a s (Mz), with twist-4 contribution, considered as 
the additional free parameters of the fit, have the following form: 

NLO HT free a s (M z ) = 0.124t{j;gjg(s*at) (35) 
±0.005{syst) ± 0.008{theory) 

NNLO HT free a s (M z ) = 0.116tg;$$(s*a*) (36) 
±0.005(syst) ± 0.004{theory) 

where the systematic uncertainties are taken from the CCFR experimental analysis, 
presented in the first work of Ref.|Q], and the theoretical uncertainties in the results 
of Eqs.(33),(35) [Eqs.(34),(36)] are estimated by the diffferences between the central 
values of the outcomes of the NNLO and NLO [N 3 LO and NNLO] fits, presented 
in Table 1, plus the arbitrariness in the application of smoothing procedure of the 
MS-scheme matching condition (which following the considerations of Ref.[Q was 
estimated as Aa s (Mz) = ±0.001) and the uncertainties in a s (Mz) due to nuclear 
effects on iron target (see discussions in Sec. 4). In estimating the theoretical errors of 
the inclusion of the N 3 LO corrections we take into account the diffferences between 
the applications of the expanded and non-expanded Pade approximants. 

It can be seen that due to the large overall number of the fitted parameters the 
results of Eqs.(35),(36) for a s (Mz) have rather large statistical uncertainties. As can 
be seen from the results of Eq.(33),(34) for the QCD coupling constant it is possible 
to decrease their values by fixing the concrete form of the twits-4 parameter h(x). 
However, if one is interested in the fitted form of the twist-4 parameter h(x), one 
should take for granted these intrinsic theoretical uncertainties of the value of a s (Mz)- 

3 (b). The results of the analysis of the experimental data: the extrac- 
tion of parameters of the twist-4 terms. 

Apart of the perturbative QCD contributions the expressions for DIS structure 
functions should contain the power-suppressed high-twist terms, which reflect the 
possible non-perturbative QCD effects. The studies of these terms have rather long 
history. At the beginning of these studies it was realized that the twist-4 contri- 
butions to structure functions should have the pole- like behavior ~ 1/((1 — x)Q 2 ) 
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p4| , [45| . This behavior was used in the phcnomcnological investigations of the earlier 
less precise than CCFR DIS vN data |4l], |47|, |48|, which together with different other 
procedures on analyzing neutrino DIS data p0| was considered as the source of 
the information about scaling violation parameters. The development of renormalon 
technique (Refs.|2(], [l9| and Ref.pH for the detailed review) pushed ahead the more 
detailed phenomenological analysis of the possibility of detecting higher-twist compo- 
nents in the available at present most precise DIS data obtained by BCDMS, SLAC, 
CCFR and other collaborations. It turned out that despite the qualitative status of 
renormalon approach the satisfactory description of the results of QCD NLO Fi SF 
analysis [jsT] in terms of IRR technique was achieved (li], ^0|. The next step was to 
clarify the status of the predictions of Ref.|ll| for the form and sign of the twist-4 
contributions to xF§ SF. The study of this problem was done in Ref.|ll| (see also 
Ref.Q). In this chapter we are discussing the results of more refined analysis of the 
behavior of the twist-4 contributions to xF^ SF at the LO, NLO, NNLO and beyond. 

In Table 3 we study the dependence of the extracted value of the parameter A 2 
from the different orders of perturbative QCD predictions, Q 2 -cuts of the CCFR ex- 
perimental data and the coefficient 7 of parton distributions model for xF$. Note, that 
the parameter A 2 was introduced in the IRR model of Eq.(28), taken from Ref . ]l9[| , 
and fixed there as A 2 ~ —0.2 GeV 2 : which is necessary for the description of the fitted 
twist-4 results of Ref. [^l| for F2 within IRR language. We found that the value of 
this parameter extracted at the LO and NLO, is negative, differ from zero for about 
one standard deviation and qualitatively in agreement with the IRR-motivated guess 
of Ref. jl9| . Moreover, the results of our LO and NLO fits are also in agreement with 
the made in Ref.|l8| earlier extraction from the old DIS neutrino data of the value 
h = —0.38 ± 0.06 GeV 2 of the different model of the twist-4 contribution to xF^, 
namely xF 3 (x, Q 2 )h/{{1 - x)Q 2 ). 

It is interesting to notice that the results of the multiloop extractions of the pa- 
rameter A 2 from the CCFR data are almost nonsensitive to the introduction of the 
additional parameter 7 in the parton distribution model. 

Another observation is that the results of Table 3 reveal that for larger Q 2 -cuts 
10 — 15 GeV 2 the results of A 2 in the LO and NLO fits are more stable to the inclusion 
of the concrete experimental data, than in the case of the low Q 2 cut (5 GeV 2 ). This 
feature can be related to the logarithmic increase of the QCD coupling constant A a at 
lower Q 2 . However, since we are interested in the extraction of the power-suppressed 
twist-4 contribution, we shall concentrate on the discussion of the more informative 
from our point of view fits with low Q 2 -c\\t 5 GeV 2 , which contain more experimental 
points and thus are more statistically motivated. 

We also observed that in the process of multiloop extraction of A 2 the tendency 
Xlo > Xnlo > Xnnlo ~ Xn3 LO takes place. It should be stressed that in opposite 
to the results of the LO and NLO extraction the NNLO value of A 2 is compatible 
with zero and is stable to the inclusion of the N 3 LO contribution to A s through the 
Pade approximants for the NS Mellin moments. 

We are now turning to the pure phenomenological extraction of the twist-4 con- 
tribution h(x) to xF% (see Eq.(24)), which is motivated by the work of Ref.j5l| for F2 
data. In the framework of this approach the x-shape of h(x) is parametrized by the 
additional parameters hi = h(xi), where xi are the points of the experimental data 
bining. The results of the multiloop extractions of these parameters are presented in 
Table 4 and are illustrated by the curves of Fig.l. 
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LO 


NLO 


NNLO 


N 3 LO 


N 3 LO(n.e.) 


X z /points 


66.2/86 


65.6/86 


65.7/86 


65.6/86 


64.3/86 


A 


5.44 ± 1.74 


3.70 ± 1.56 


4.54 ± 0.88 


5.24 ± 0.92 


5.16 ± 0.75 


b 


0.74 ± 0.10 


0.66 ± 0.11 


0.69 ± 0.06 


0.72 ± 0.06 


0.73± 0.05 


c 


4.00 ± 0.18 


3.78 ± 0.21 


3.72 ± 0.19 


3.58 ± 0.26 


3.48 ± 0.26 


7 


1.72 ± 1.25 


2.86 ± 1.72 


1.43 ± 0.69 


0.72 ± 0.63 


0.71 ± 0.51 


A/T C 
1V1 J 

[MeV] 


338 ± 169 


428 ± 158 


264 ± 85 


248 ± 76 


310 ± 100 














h{x 4 ) [GeV 2 ] 


0.0125 


0.206 ± 0.321 


0.213 ± 0.332 


0.170 ± 0.302 


0.165 ± 0.300 


0.190 ± 0.303 


0.0175 


0.061 ± 0.268 


0.091 ± 0.289 


0.030 ± 0.241 


0.020 ± 0.240 


0.051 ± 0.247 


0.025 


0.146 ± 0.204 


0.220 ± 0.241 


0.136 ± 0.178 


0.118 ± 0.176 


0.156 ± 0.186 


0.035 


-0.021 ± 0.185 


0.114 ± 0.240 


0.008 ± 0.179 


-0.021 ± 0.175 


0.023 ± 0.182 


0.050 


0.031 ± 0.142 


0.245 ± 0.230 


0.124 ± 0.164 


0.080 ± 0.154 


0.130 ± 0.162 


0.070 


-0.145 ± 0.127 


0.138 ± 0.233 


0.033 ± 0.165 


-0.026 ± 0.146 


0.028 ± 0.163 


0.090 


-0.177 ± 0.125 


0.139 ± 0.225 


0.076 ± 0.166 


0.008 ± 0.143 


0.063 ± 0.161 


0.110 


-0.340 ± 0.126 


-0.015 ± 0.205 


-0.026 ± 0.166 


-0.095 ± 0.145 


-0.042 ± 0.163 


0.140 


-0.404 ±0.114 


-0.092 ± 0.147 


-0.027 ± 0.141 


-0.088 ± 0.128 


-0.048 ± 0.137 


0.180 


-0.350 ± 0.164 


-0.077 ± 0.122 


0.054 ± 0.127 


0.016 ± 0.130 


0.021 ± 0.128 


0.225 


-0.554 ± 0.237 


-0.348 ± 0.170 


-0.167 ± 0.135 


-0.174 ± 0.139 


-0.221 ± 0.142 


0.275 


-0.563 ± 0.334 


-0.462 ± 0.272 


-0.196 ± 0.193 


-0.171 ± 0.189 


-0.264 ± 0.206 


0.350 


-0.314 ± 0.418 


-0.368 ± 0.371 


0.070 ± 0.200 


0.128 ± 0.188 


0.026 ± 0.206 


0.450 


-0.117 ± 0.415 


-0.266 ± 0.401 


0.183 ± 0.213 


0.237 ± 0.191 


0.119 ± 0.209 


0.550 


0.087 ± 0.333 


-0.109 ± 0.352 


0.097 ± 0.234 


0.116 ± 0.224 


-0.051 ± 0.278 


0.650 


0.377 ± 0.215 


0.221 ± 0.244 


0.259 ± 0.188 


0.271 ± 0.189 


0.184 ± 0.220 



Table 4. The results of the LO, NLO (N max = 10), NNLO and the N 3 LO {N Max = 6) expanded 
and non-expanded (n.e.) Pade QCD fit (with TMC) of the CCFR'97 xF 3 SF data for the values of 
HT contributions h(x) and for the parameters A, 6, c, 7 with the corresponding statistical errors. 



Looking carefully on Table 4 and Fig.l we observe the following features: 

1. The x-shape of the twist-4 parameter is not inconsistent with the expected rise 
of h(x) for x — > 1 1 44, 45) in all orders of perturbation theory; 

2. The values of the parameters h{x{) at the upper and lower points of kinematic 
region (aii6=0.650 and xi=0.0125) are stable to the inclusion of the higher order 
perturbative QCD corrections and application of the Pade resummation tech- 
nique. At large values of x this feature is in agreement with the previous state- 
ment; 

3. The function h(x) seems to cross zero twice: at small x of order 0.03 and larger 
x about 0.4. It should be noted that the sign-alternating behavior of the twist-4 
contributions to DIS structure functions was qualitatively predicted in Ref . |47j ] ; 

4. In the LO and NLO our results are in qualitative agreement with the IRR pre- 
diction of Ref. |l9| (for discussions see Ref. ^lj ) ; 

5. In the NNLO this agreement is not so obvious, though the certain tendency of 
following the general shape of IRR prediction llll still survives; 
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Fig.l The results of the LO, NLO, NNLO and N 3 LO [0/2] Pade extarctions of 
the twist-4 contributions of h(x). The solid line is the IRR- model prediction of 

Ref.|§. 

6. However, at the NNLO we observe the partial nullification of h{x) within statis- 
tics error bars. Thus we conclude that the inclusion of the NNLO corrections 
into the game is shadowing the effects of the power suppressed terms, or that the 
effects of the twist-4 corrections are nondetectable at the NNLO. This property 
was previously observed at the LO as the result of the less precise DIS neutrino 
data in Ref. fill . In the modern experimental situation, namely in the process of 
the analysis of the more precise DIS neutrino data of CCFR collaboration, we 
are observing this feature at the NNLO; 

7. We checked the reliability of the foundation of partial nullification of h(x) at the 
NNLO by going beyond this perturbative approximation using the methods of 
Pade approximants. The result of this analysis reveals the stability of the NNLO 
results for h(x) and its partial nullification; 

8. The property of partial nullification of h(x) at the NNLO and N 3 LO is identical 
to the effect of nullification of the IRR model parameter A 2 at the NNLO and 
N 3 LO (see Table 3); 

9. These observed properties clarifies why the results of the NNLO and N 3 LO fits 
for A-^L. , presented in Tables 3, practically do not depend from the inclusion of 
the twist-4 contribution through the IRR model. Indeed, at this level the twist-4 
terms have almost zero effect. 
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To our point of view the foundations (8)- (9) reflects the selfconsistency of the 
results of our different fits with twist-4 included by different ways. 

4. The quest of the inclusion of the effects of nuclear corrections. 

The effects of nuclear corrections are remaining the important source of the uncer- 
tainties of the analysis of the D1S data. This is especially important for the experi- 
ments on heavy targets and in the case of CCFR data-on iron 56 Fe. 

The attempts to study these effects were done in Ref. |52] in the framework of 
Deutron-motivated model. The satisfactory QCD description of the CCFR data for 
xF s was achieved due to the reason that in this case the nuclear effects do not exceed 5 
% effect. However, the more realistic description of nuclear effect for xF 3 on 56 Fe J53| 
revealed the appearance of new 1/Q 2 and 1/M corrections for NS moments (where M 
is the mass of the nucleon), which have the following form 

m£(q 2 )/a = ^i + ^(„-i) + ±E!z ? ^-i) + o(-L^M^(g 2 ) 

+ < Ap 2 > d p 2M» s (Q 2 ) 



2 < > 
3Q 2 



n(n+l)M^ 2 (Q 2 ) (37) 



where for 56 Fe the parameters of the nuclear model, adopted in Ref. |53|| are < e >« 
-56 MeV, < p 2 > /(2M) » 35 MeV, < Ap 2 > Fe « -0.17 GeV and the derivative 
d p 2M n (Q 2 ) is taking into account that the target momentum p can be generally of- 
mass-shell. This effect is resulting in the following independent from the nuclear 
content contribution JIH 

d p iM n (Q 2 ) = dpM™ + J (m^ S + M 2 d p2 M^ (38) 

where the numerical values of t) p 2M" were also presented in Ref. [[53f . 

Note, that the effects of the nuclear corrections in DIS were also recently studied 
in Ref. f| in the case of xF 3 SF and in Ref.||,|6) in the case of F 2 SF (for the 
earlier related works see e.g. Ref.j57j]). However, in our studies we will concentrate 
ourselves on the consideration of the results of Ref. |53| . 

We included the corrections of Eqs.(37)-(38) into our fits and observed the unac- 
ceptable increase of \ 2 value. We think that this can be related to to the manifestation 
of the possible asymptotic character of the 1 /M-expansion in Eq. (37) , since the third 
term in the brackets of the r.h.s. of Eq.(32) becomes comparable with the first term 
(which is equal to unit) for the n ~ 8 used in our fits. Note that the moments 
with large n are important in the reconstruction of the behavior of the xF 3 SF at 
x — ► 1. This observed feature necessitates the derivation of the explicit expression for 
M^(Q 2 ), which is not expanded in powers of l/M-terms. It should be noted that the 
problem of the possible asymptotic nature of the power suppressed expansions was 
mentioned in the case of Ellis- Jaffe and Bjorken DIS sum rules in Ref. pSj] . 

Another possibility of the nonconvergence of our fits with nuclear corrections of 
Eq.(37) taken into account might be related to the fact that the parton distribution 
model for the nuclear SF xF 3 Fe can be different from the canonical model, used by 
us. 

In any case we think that the problem of the taking into account of the heavy 
nuclear effects in the process of fits of xF 3 data is really on the agenda. 
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Conclusion 

In this work we presented the results of the extractions of a s (Mz) and twist-4 terms 
from the QCD analysis of the CCFR data with taking into account definite QCD 
corrections at the NNLO and beyond. Within experimental and theoretical errors 
our results for a s (Mz) are in agreement with other extractions of this fundamental 
parameter, including its world average value a s (Mz) = 0.118 ± 0.005. 

Our estimate of the NNLO theoretical uncertainties is based on application of the 
[0/2] Pade approach at the N 3 LO level. The uncertainties of our NNLO analysis can 
be decreased after explicit NNLO calculations of the NS Altarelli-Parisi kernel. 

As to the twist-4 terms, we found that despite the qualitative agreement with the 
IRR model prediction, at the NNLO level they have the tendency to decrease and are 
stable to the application of the [0/2] Pade motivated N 3 LO analysis. 

This feature can be related to the fact that the analysis of the CCFR data can not 
distinguish the twist-4 1/Q 2 terms from the NNLO perurbative QCD approximations 
of the Mellin moments. This possible explanation lies in the lines of the results of the 
LO analysis of the old less precise neutrino DIS data, made by the authors of Ref. p§] , 
who were unable to distinguish between LO logarithmic and l/Q 2 -behavior of the 
QCD contributions to Mellin moments of xF 3 . The achieved in our days increase of 
experimental precision might move this effect to the NNLO. 

Another related explanation is that the observed by us NNLO effect is manifesting 
itself in view of the fact that the detected by us twist-4 terms come from the partial 
summation of the definite terms of the asymptotic perturbative QCD series and thus 
the increase of the order of perturbative QCD analysis effectively suppresses the re- 
maining sum of the perturbative QCD contribution. Unfortunately, at the NNLO level 
we can not detect the true twist-4 terms, which reflect the non-trivial non-perturbative 
nature of QCD vacuum. One can hope that the future experiments of NuTeV collab- 
oration will allow to get the new experimental data at the precision level, necessary 
for extracting more detailed information about higher twist contributions to structure 
functions and will help to clarify the reason of the disagreement of the low x CCFR 
data for F 2 SF with the ones, obtained by the BCDMS collaboration. We hope to 
return to the NNLO analysis of the experimental data of BCDMS collaboration in the 
nearest future. 

Acknowledgments 

The part of this work was done when one of us (AVS) was visiting Santiago de 
Compostela Univ. He is grateful to his colleagues for hospitality in Spain. Our studies 
were continued when ALK was participating at the LHC meeting at CERN (March, 
1998). He is grateful to the Organizers of this Meeting for the invitation and support. 

We are grateful to A.V.Kotikov for participation at the first stage of these our 
studies and for discussions. The useful comments of S.I. Alekhin are gratefully ac- 
knowledged. Special thanks are due to S.A.Kulagin for sharing his points of view on 
the currents status of the studies of the nuclear effects in DIS. 

The preliminary versions of this work were reported at Quarks-98 Int. Seminar 
(Suzdal, May, 1998), Small- a; physics and Light Front Dynamics in QCD Conference 
(St. Petersburg, July, 1998) and XI Int. Conference on Problems of Quantum Field 
Theory in memory of D.I. Blokhintsev (Dubna, July, 1998). We would to thank 
the participants of these conferences, and especially S.J. Brodsky, V.M. Braun, K.G. 
Chetyrkin, B.L. Ioffe, L.N. Lipatov, E.A. Kuraev and A.V. Radyushkin for the interest 
in this work. 



18 



This work is supported by the Russian Fund for Fundamental Research, Grant N 
96-02-18897. The work of G.P. was supported by CICYT (Grant N AEN96-1773) and 
Xunta de Galicia (Grant N XUGA-20604A96). 



References 

CCFRR Collab., D. MacFarlane et a!., Z. Phys. C26 (1984) 1. 

W.G. Seligman, Columbia Univ, Thesis R-1257, CU-368, Nevis- 292, 1997; 
CCFR-NuTeV Collab., W.G. Seligman et ah, Phys. Rev. Lett. 79 (1997) 1213. 

G. Altarelli and G. Parisi, Nucl. Phys. B126 (1977) 298; V.N. Gribov and L.N. Li- 
patov, Sov. J. Nucl. Phys. 15 (1972) 438; L.N. Lipatov, Sov. J. Nucl. Phys. 20 
(1975) 94; Yu.L. Dokshitzer, JETP 46 (1977) 641. 

S. Kuhlmann, HI . Lai and W.K. Tung, preprint MSU-HEP-70316; CTEQ-705 
( |hcp-ph/9704338| ) 



A.D. Martin, R.C. Roberts, W.J. Stirling and R.S. Thorne, Eur. Phys. J. C4 
(1998) 463. 



M. Gluck, E .Reya and A. Vogt, preprint DO-TH 98/07; WUE-ITP-98-019 flhep- 
ph/9806404| ). 



A.L. Kataev and A.V. Sidorov, Phys. Lett. B331 (1994) 179. 
J. Chyla and J. Rames, Phys. Lett. B343 (1995) 351. 

A.L. Kataev, A.V. Kotikov, G. Parente and A.V. Sidorov, Phys. Lett. B388 
(1996) 179. 

A.V. Sidorov, Phys. Lett. B389 (1996) 379; 

A.V. Sidorov, JINR Rapid Commun. 80 (1996) 11. 

A.L. Kataev, A.V. Kotikov, G. Parente and A.V. Sidorov, Phys. Lett. B417 
(1998) 374. 

A.L. Kataev, A.V. Kotikov, G. Parente and A.V. Sidorov, Proc. QCD-97 Work- 



shop, Montp ellier, July 1997; Nucl. Phys. B Proc. Suppl. 64 (1998) 138 flhep- 
ph/9709509| ). 



G. Parisi and N. Sourlas, Nucl. Phys. B151 (1979) 421. 
J. Chyla and J. Ramez, Z. Phys. C31 (1986) 151. 

V.G. Krivokhizhin, S.P. Kurlovich, V.V. Sanadze, LA. Savin, A.V. Sidorov and 
N.B. Skachkov Z. Phys. C36 (1987) 51; 

V.G. Krivokhizhin, S.P. Kurlovich, R. Lednicky, S. Nemechek, V.V. Sanadze, LA. 
Savin, A.V. Sidorov and N.B. Skachkov, Z. Phys. C48 (1990) 347. 

E.B. Zijlstra and W.L. van Neerven, Nucl. Phys. B417 (1994) 61. 

S.A. Larin, T. van Ritbergen and J. A.M. Vermaseren, Nucl. Phys. B427 (1994) 
41; 

S.A. Larin, P. Nogueira, T. van Ritbergen and J. A.M. Vermaseren, Nucl. Phys. 
B492 (1997) 338. 

G. Parente, A.V. Kotikov and V.G. Krivokhizhin, Phys. Lett. B333 (1994) 190. 
M. Dasgupta and B.R. Webber, Phys. Lett. B333 (1994) 190. 



19 



[20] M. Beneke and V.M. Braun, Nucl. Phys. B454 (1995) 253; 

Yu.L. Dokshitzer, G. Marchesini and B.R. Webber, Nu cl. Phys. B469 ( 1996) 93; 
R. Akhoury and V.I. Zakharov, report UM-TH-97-05 ( |hep-ph/9701378[ ); 
R. Akhoury and V.I. Zakharov, Proc, QCD- 97 Workshop, Mo ntpellier, July 1998, 
Nucl. Phys. B. Proc. Suppl. 64 (1998) 350 ( |hep-ph/9710257D ; 
M. Maul, E. Stein, L. Mankiewicz, M. Mayer-Hermann and A. Schaier, Invited 
talk at 2nd Topical Workshop of DIS of pola rized Targets: T heory Meets Exper- 
iment, Zcuthcn, Germany, September 1997 ( hcp-ph/9710392 ). 

M. Beneke, preprint CERN-TH-98-233 flhep-ph/9807443| ). 

W. Bernreuther and W. Wetzel, Nucl. Phys. B197 (1982) 228. 

W. Marciano, Phys. Rev. D29 (1984) 580. 

D.V. Shirkov and S.V. Mikhailov, Z. Phys. C63 (1994) 463; 
D.V. Shirkov, A.V. Sidorov and S.V. Mikhailov, jhep-ph/9707524 . 

S.A. Larin, T. van Ritbergen and J. A.M. Vermaseren, Nucl. Phys. B438 (1995) 
278; for the definite cross-checks see M. Spira, Fortsch. Phys. 46 (1998) 203. 

T. van Ritbergen, J. A.M. Vermaseren and S.A. Larin, Phys. Lett. B400 (1997) 
397. 

K.G. Chetyrkin, B.A. Kniehl and M. Steinhauser, Phys. Rev. Lett. 74 (1995) 
2184. 

G. A. Baker, Jr. and P. Graves-Morris, Pade approximants, Addison- Wesley Pub. 
Co., 1981; 

CM. Bender and S.A. Orszag, Advanced mathematical methods for scientists 
and engineers, McGraw-Hill Book Company, 1978. 

MA. Samuel, J. Ellis and M. Karliner, Phys. Rev.Lett. 74 (1995) 4380. 

J. Ellis, E. Gardi, M. Karliner and M.A. Samuel, Phys. Lett. B366 (1996) 268. 

A.V. Sidorov, Proc. Int. Conf. Hadron Structure-91, Stara Lesna, Czecho- 
slovakia, 1991, Physics and Applications, Kosice, 1992, p. 230. 

A.L. Kataev and V.V. Starshenko, Mod. Phys. Lett. A10 (1995) 235. 

J. Ellis, M. Karliner and M.A. Samuel, Phys. Lett. B400 (1997) 176. 

J. Ellis, I. Jack, D.R.T. Jones, M. Karliner and M.A. Samuel, Phys. Rev. D57 
(1998) 2665; 

V. Elias, T.G. Steele, F. Chishtie, R. Migneron and K. Sprague, jhcp-ph/9806324 . 
K.G. Chetyrkin, Phys. Lett. B404 (1997) 161; 

J. A. Vermaseren, S.A. Larin and T. van Ritbergen, Phys. Lett. B405 (1997) 327. 

CM. Bender and S. Bocttcher, J. Math. Phys. 35 (1994) 1914; 
CM. Bender, S. Boettcher and L.N. Lipatov, Phys. Rev. Lett. 65 (1992) 3674; 
Phys. Rev. D46 (1992) 5557. 

A.A.Pcnin and A. A. Pivovarov, Phys. Lett. B357 (1995) 427. 
L. Durand and G. Jaczko, hcp-ph/9804362 . 

H. Georgi and H.D. Politzer, Phys. Rev. D14 (1976) 1829. 

S.J. Brodsky, J. Ellis, E. Gardi, M. Karliner and M.A. Samuel, Phys. Rev. D56 
(1997) 6980; 

G. Cvetic, Nucl. Phys. B517 (1998) 506; 
G. Cvetic and R. Kogerler, Nucl. Phys. B522 (1998) 396. 



20 



J. Chyla and A.L.Kataev, Phys.Lett. B297 (1992) 385. 

J. Chyla and A.L.Kataev, in Report ol Working Group on Precisison calculations 
for the Z°-resonance, Eds. D.Yu. Bardin, W. Hollik and G.Passarino; CERN 
Yellow Report 95-03, p.313 Qhcp-ph/9502383| ). 

A. V. Kotikov, G. Parente and J. Sanchez Guillen, Z. Phys. C58 (1993) 465. 
E.Berger and S.J. Brodsky, Phys. Rev. Lett 42 (1979) 940. 

J. Gunion, P. Nason and R. Banknbecler, Phys. Rev. D29 (1984) 2491. 

L.F. Abbot and R.M. Barnett, Ann. of Phys. 125 (1980) 276. 

V.A. Bednyakov, P.S. Isaev, S.G. Kovalenko, Yad. Fiz. 40 (1984) 770 (Sov. J. 
Nucl. Phys. 40 (1984) 494). 

B. R. Martin and G. Shaw, Z. Phys. C33 (1986) 99-106. 

A. De Rujula, H. Georgi and H.D. Poitzer, Phys. Rev D10 (1974) 2141. 

A. Gonzales-Arroya, C. Lopez and F.J. Yndurain, Nucl. Phys. B153 (1979) 161. 

M. Virchaux and A. Milsztajn, Phys. Lett. B274 (1992) 221. 

A.V. Sidorov and M.V. Tokarev, Phys. Lett. 358 (1995) 353; 
M.V. Tokarev and A.V. Sidorov, II Nuovo Cim. A110 (1997) 1401. 

S.A. Kulagin, preprint INR P0967/97 ( |nucl-th/9801035[ ), to be published in Nucl. 
Phys. A 



E. Marco E. Oset and S.K. Singh, preprint FTUV-98-35 (1998) ( |nucl- 
th/9804068D. 



K.J. Eskola V.J. Kolhinen and C.A. Salgado, preprint JYFL-8-38 (1998) flhep- 
ph/9807297l ). 



CD. Cothran, D. B. Day and S. Liuti, Phys. Lett. B421 (1998) 46. 

SV. Akulinichev, S.A. Kulagin and G.M. Vagradov, Phys. Lett. B158 (1985) 
485; 

SV. Akulinichev, S. Shlomo, S.A. Kulagin and G.M. Vagradov, Phys. Rev. Lett. 
55 (1985) 2239; 

B.I. Birbrair, EM. Levin and A.V. Shuvaev, Nucl. Phys. A491 (1985) 619. 

B.L. Ioffe, L ectures at Int. S chool of Nucleon Structure, Erice, 1995, Preprint 
ITEP-62-95 ( |hep-ph/951140ll ). 

S.A. Kulagin, private communication, September 1998. 



21 



